Multiple quantum NMR dynamics of spin-^ carrying molecules of a gas in nanopores 
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We consider the multiple quantum (MQ) NMR dynamics of a gas of spin carrying molecules in 
nanocavities. MQ NMR dynamics is determined by the residual dipole-dipole interactions which 
are not averaged completely due to the molecular diffusion in nanopores. Since the averaged non- 
secular Hamiltonian describing MQ NMR dynamics depends on only one coupling constant, this 
Hamiltonian commutes with the square of the total spin angular momentum I 2 . We use the basis of 
common eigenstates of I 2 and the projection of / on the external magnetic field for investigation of 
MQ NMR dynamics. This approach allows us to study MQ NMR dynamics in systems consisting of 
several hundreds of spins. The analytical approximation of the stationary profile of MQ coherences 
is obtained. The analytical expressions for MQ NMR coherence intensities of the five-spin system 
in a nanopore are found. Numerical investigations allow us to find the dependencies of intensities 
of MQ coherences on their orders (the profiles of MQ coherences) in systems consisting of 600 spins 
and even more. It is shown that the stationary MQ coherence profile in the considered system is an 
exponential one. 



I. INTRODUCTION 

Multiple quantum (MQ) NMR dynamics is a basis of 
MQ NMR spectroscopy [l| which, in turn, is a powerful 
tool to study the nuclear spin distributions in different 
systems (for example, in liquid crystals [2j], simple or- 
ganic systems [1], amorphous hydrogenated silicon [3[, 
etc.). MQ NMR has been used [4|, [5[ to investigate the 
size of spin clusters when the growth of MQ clusters oc- 
curs in the process of the irradiation of the spin system on 
the preparation period of the MQ NMR experiment 
The unique possibilities of MQ NMR to study dynamics 
of many-spin clusters have been recently used 0, Q for 
a measurement of the decoherence rate for highly cor- 
related spin states. The scaling of the decoherence rate 
with the number of correlated spins has been also pre- 
sented 0]. 

The theoretical description of MQ NMR dynamics is 
a very difficult task because this problem is a many- 
spin and multiple-quantum one. Although the first phe- 
nomenological approach was developed together with the 
experimental realization of MQ NMR [1] the consistent 
quantum-mechanical theory is not developed up to now 
except for one-dimensional systems. In particular, the 
one-dimensional spin chain with nearest neighbor double 
quantum Hamiltonian [1] is exactly solvable and it has 
been shown that, starting with a thermodynamic equilib- 
rium state, only zero and double quantum coherences are 
produced @-[l0]. This conclusion has been confirmed ex- 
perimentally for relatively small excitation times [ll|, [l2[ • 
However, next-nearest couplings and other distant inter- 
actions lead to higher order coherences in MQ NMR spec- 
tra. These interactions are beyond the exact solvable 
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models [l|-[lo|- Thus one should apply numerical meth- 
ods in order to take into account distant spin-spin cou- 
plings. However, even supercomputer calculations allow 
us to study MQ NMR dynamics of a spin chain consisting 
of no more than fifteen spins [l3[ , which is insufficient to 
solve some subtle problems of MQ NMR dynamics. In 
particular, the dependence of the MQ coherence inten- 
sities on their orders (the profile of MQ coherences [J]) 
can't be found. 

Nanosize systems exhibit new possibilities for investi- 
gations of MQ NMR dynamics. It is well known [14] that 
dipole-dipole interactions (DDI) of spin carrying atoms 
(molecules) of a gas in the non-spherical nanopores in 
a strong external magnetic field are not averaged com- 
pletely due to the molecular diffusion [l4], EH . We em- 
phasize that the residual averaged DDI are determined 
by only one coupling constant which is the same for all 
pairs of interacting spins [H], EH- As a resu lt, the av- 
eraged non-secular two-spin/ two-quantum Hamiltonian 
[l| describing MQ NMR dynamics commutes with the 
operator of the square of the total spin angular momen- 
tum I 2 . Note that eigenstates of the projection of the 
total spin angular momentum on the external magnetic 
field I z are usually used as a basis (the multiplicative 
basis) in order to describe MQ NMR dynamics. Since 
[Pjz] = 0, it is suitable to study MQ NMR dynamics 
in a nanopore in the basis of common eigenstates of I 2 
and I z . Since there are no transitions with changing I 2 in 
the MQ NMR experiments the problem splits into a set of 
simpler problems for different values of I 2 and it is possi- 
ble to solve the problem of the exponential growth of the 
Hilbert space dimension with an increase of a number of 
spins. For this reason one can investigate MQ NMR dy- 
namics in systems consisting of several hundreds of spins 
and solve the problem of the profile of intensities of MQ 
NMR coherences. 

The present paper is devoted to the investigation of 
MQ NMR dynamics of spins coupled by the DDI in a 
nanopore. The paper is organized as follows. In Sec. II, 
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spin- 1/2 probe 




FIG. 1: Spin-carrying molecules (atoms) in a nanopore in an 
external magnetic field. 



the theory of MQ NMR dynamics of spins coupled by 
the DDI is developed. The numerical algorithm for MQ 
NMR dynamics in a nanopore is also discussed in Sec. 
II. The analytical solution for MQ NMR dynamics of a 
five-spin system is obtained in Sec. III. The profiles of 
intensities of MQ NMR coherences for the systems con- 
sisting of 200 -T- 600 spins are presented in Sec. IV. Here 
we give also analytical approximations of the intensity 
profiles with smooth curves and compare them with the 
simple physical estimations. We briefly summarize our re- 
sults and discuss further perspectives in the concluding 
Sec. V. 



II. MQ NMR DYNAMICS OF SPINS IN A 
NANOPORE 

We consider a system of nuclear spins (s — 1/2) cou- 
pled by the DDI in a strong external magnetic field in a 
nanopore [Fig. [T]. The secular part of the DDI Hamilto- 
nian [16[ has the following form: 



Djk[2ij Z ikz 

j<k 



(1) 



where D 



jk 



2r 



jk 



(1 — 3 cos 2 Ojk) is the coupling constant 

between spins j and /c, 7 is the gyromagnetic ratio, rjk 
is the distance between spins j and /c, and Ojk is the an- 
gle between the internuclear vector fjk and the external 
magnetic field which is directed along z axis. The op- 
erator Ij a (a = x, y, z) is the projection of the angular 
spin momentum operator on the axis a; 1^ and IJ are 
the raising and lowering operators of spin j. 

The standard MQ NMR experiment consists of four 
distinct periods of time [Fig. [2]: preparation (r), evolu- 



tion (£1), mixing (r) and detection (t 2 ) Q- MQ coher- 
ences are created by the multipulse sequence consisting of 
eight-pulse cycles on the preparation period [l|. In the 
rotating reference frame [16], the averaged non-secular 
two-spin /two-quantum Hamiltonian, Hmq, describing 
MQ dynamics at the preparation period can be written 
as 



MQ 



ff(2) +fr (-2) 



where 



(2) 



(3) 



j<k 



Since the time of the molecular diffusion in nanopores 
is much shorter than both the dipolar time t ~ uj^ c 
(uof oc = Tr(Hj z /Tr(I%))^ and the period of the multi- 
pulse sequence at the preparation period of the MQ NMR 
experiment [l| , one can assume that the spin dynamics is 
governed by the averaged dipolar coupling constant, D, 
which is the same for all spin pairs. Although the prob- 
lem is described by only one averaged coupling constant 
it is the many-spin one with the Hamiltonian 

H MQ = -j{(I + ) 2 + (I-) 2 }, (4) 

N 

^ — ^2 it an d N is the number of spins in the 



where / 



nanopore. 

In order to investigate MQ NMR dynamics of the sys- 
tem one should find the density matrix p solving the Li- 
ouville evolution equation [16[ (h=l): 



.dp 
dr 



(5) 



with the initial density matrix p(0) = I z in the high 
temperature approximation [16[. It proves convenient to 
expand the spin density matrix, p, in the series as 



(6) 



where pkij) is the contribution to p(r) from MQ coher- 
ences of the kth order. Then the intensity Jfe(r) of the 
MQ coherence of order A:, which is an observable in the 
MQ NMR experiment [l|, is determined as 



Jfe(r) 



Tr{/? fc (r) • p-fe(r)} 
Tr(/f) 



(7) 



The Hamiltonian Hmq leads to the emergence of MQ 
coherences of even orders only, and the coherence order 
cannot exceed the number of spins N [l|. 

Note that the eigenvalues and the eigenvectors of Hmq 
have the following property. Let A and u be the eigen- 
value and the corresponding eigenvector of the Hamilto- 
nian Hmq- Since 



(8) 
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~ H MQ 
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FIG. 2: The basic scheme of the MQ NMR experiment. 



one can find that —A and exp(— i^I z )u are also the eigen- 
value and the eigenvector of Hmq- 

Since the square of the total spin angular momentum 
I 2 commutes with projections of / on an arbitrary direc- 
tion, we have from Eq. (j4j) that 



[H MQ J 2 ]=0. 



(9) 



Hereafter we will use the basis consisting of the common 
eigenstates I 2 and I z to study MQ NMR dynamics. In 
this basis, the Hamiltonian Hmq consists of the different 

5 

blocks H M q corresponding to the different values of the 

total spin angular momentum S (S 2 = S(S + 1), S = 
— 1,^ — 2,...,^ — [^r] ; [i] is an integer part of i): 



N N 
2 ' 2 



—S + 1, —S + 2, . . . , S — 1, S. Hence one can 

{M\(I + ) 2 \M - 2) = (M - 2\(r) 2 \M) (15) 
= y/(S + M)(S + M — 1)(S - M + 1)(5 - M + 2), 

where M = -5 + 2, -5 + 3, . . . , S - 1, 5 and all other 

elements are zero. 

Eq. (fT5|) determines the matrix representation of the 
g 

block H M q. The dimension of this block is 25 + 1. The 
total dimension of the Hamiltonian, 2^, is determined by 
the sum of the dimensions of all the blocks as follows 



^njv(5)(25+l) 



N\(2S+ l) 2 



+ 5 + -5)! 



(16) 



The proof of the relationship (fTB|) is given in Appendix. 

Now we use one more symmetry of Hmq- Since the 
Hamiltonian Hmq and all its blocks are invariant with 
respect to the rotation by the angle 7r about the z axis, 
the operator exp(z7r/ 2 ) is an integral of motion and 



H 



MQ 



dmg{H MQ ,H MQ 



MQ 



J }. (10) 



\H M Q,exp(inIz)] = 0. 



(17) 



Taking into account the block structure of the Hamil- 
tonian Hmq and the diagonal structure of I z one can 
find that the density matrix p{r) consists of blocks p s (r) 
(S = f , f - 1, . . . , f - [f ] ) as well. We will denote as 
pf (t) the contribution to p S (r) from MQ coherence of 
order k. Then the contribution Jk,s(j) to the intensity 
of the /c-th order MQ NMR coherence is determined as 



Jk,s(r) 



P S - k ) 



Tr(/f) 



(11) 



It means that the problem is reduced to the set of the 

analogous tasks for each block H M q . The number of the 
energy levels, un(S), for the total momentum S in the 
TV-spin system is [17| 



It means that the 2^ x 2^ Hamiltonian matrix is reduced 

to two 2 Ar_1 x 2 Ar_1 submatrices. The same is valid for all 
g 

blocks H M q which split into two subblocks: H M q and 

5 

H M q. This reduction is valid for arbitrary N. However 
for odd TV, both submatrices yield the same contribution 
into the profile of the MQ NMR coherences [13[. Con- 
sequently, one should solve the problem using only one 



->iV-l 



x 2 



N-l 



submatrix and double the obtained inten- 



sities [13]. Below all calculations are performed when 
the number of spins is odd. The numerical algorithm for 
NMR spin dynamics with the given total spin angular 
momentum is based on the diagonalization of the blocks 

Hmq (0 < *S < ^ of the Hamiltonian Hmq and is 
described in Ref. KLOL 



n N (S) 



N\(2S + 1) 



(f + 5 + l)!(f -5)!' 



0<S'< 



N 



(12) 



which is also the multiplicity of the intensities Jk,s( r )- 
Then the observable intensities of MQ NMR coherences 
J k (r) (—TV <k<N) are 



J k(r) = ^2n N (S)J ky s(r). 



(13) 



In order to find expressions for J k (r) we need the matrix 
representations of J ± . Non-zero elements of these repre- 
sentations related with the total spin angular momentum 
S are [T3| 

(M|J+|M-1) = (\M-l\r\M) 



= V(5 + M)(5- M + 1), (14) 



III. THE EXACT SOLUTION FOR MQ NMR 
DYNAMICS OF A FIVE-SPIN SYSTEM IN A 
NANOPORE. 

We consider a system of N = 5 spins coupled by the 

averaged DDI in a nanopore. The values S of the total 

spin angular momentum are 5/2, 3/2 and 1/2 for this 

case. It was mentioned in Sec. II that if the number of 

spins is odd then each block H M q splits into two ones 
§ 

H M q and H M q and we can consider only one of them in 
order to calculate the intensities of the MQ NMR coher- 
ences [l3|| . Using Eq. ([T5|) one can obtain the eigenvalues 



l_5 y2 

1,2,3) of H M q which are the following 
= -y/7D, A$ 2 = V7D, 



\ (!) 
A 5/2 



x(3) 
A 5/2 



0. 



(18) 
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The appropriate set of eigenvectors reads as follows 



u 



u 



(1) = 

5/2 

(2) _ 
5/2 — 

(3) = 
5/2 



1 5 1 3_ 

1/5 J 3_\ 

2V 7' V2'2V7) 

,0, 



(19) 



VT4 



Similarly, the eigenvalues and eigenvectors of the block 
H~mq 2 are following 



X 3/2 
X 3/2 



4a 



A 



(2) 



2 " > ■ "3/2 

- -V 



2 ' 

(2) 
^3/2 — 



(20) 



1 1 

71' 71 



|-l/2 

The block H M q is a scalar, 



T7+V2 



MQ 



0. 



(21) 



The solutions of Eq. ([5]), pn 5 ^ 2 (^ = 1,3,5), where the 

Hamiltonian Hmq is changed by the Hamiltonian H M q , 
is 



P +n/2 (r) = U +n/2 e 



iA 



+ »/2„ 



U +n/2 



+n/2 JT JA+ n / 2 T n + 



+n/2> 



(22) 



where A +n / 2 is the diagonal matrix of eigenvalues and 



F +n/2 



U +n /2 is the matrix of eigenvectors of the block H M q 

(n = 1,3,5), and the initial density matrices p^ n ^ 2 are 
the following 



where 
an 

«22 
«23 
«33 



98 L 

i /To 
rVT 

24 
~49 
1 

14 



36 cos(V7Dr) - 2 cos(2 VlDr) + 15 
^ cos (V^Dr) -9 sin(v / 7^r), 



— \/5 sin 4 



V7Di 



4cos(2v / 7Dr) +3 



(25) 



3 2 

~7 V 7 
3 

98 



z cos (V^Dr) +5 sin(v / 7Dr), 
60cos(v / 7^r) + 6 cos(2 V7i>r) - 17 



The analogous calculations for the matrix p +3 / 2 using 
Eqs. ([20I22123D with n = 3 yield 



cos ( V^r) + \ -i sin( v / 3^r) 
z sin 



.(26) 



o+ 3 / 2 



Finally, Eq.([22]) with n = 1 in view of Eqs. (ZD and (|23l 
yields 



P +1/2 - Po +1/2 = ± 



(27) 



Only MQ coherences of zeroth, plus/minus second, and 
plus/minus fourth orders appear in the considered sys- 
tem. These intensities can be calculated with Eqs. (0, 

(EUEZJ): 



Mr) 



20 



i=l 



L(„ 5 (5/2)5>f | 2 (28) 

^s(l/2)) = 
27825 + 9604 • cos(2\/3£»t) 



+n 5 (3/2)5>: 

i=l 

1 



3/2.2 



48020 

+ 2520cos(v / 7Dr) + 7560 cos(2v / 7^r) 
+ 360 cos(3v / 7^r) + 151 cos(4v / 7^ ) r) 



p +1/2 = l/2, p +3/2 =( 3 / 2 ^3), (23) J ±2 (r) 
5/2 

+5/2 =[01/2 







-3/2 



After calculations using Eqs. (|18I19I22I23|) with n = 5 
one obtains 



P +5/2 (r) 



an za i2 a 13 
-ia 12 a 2 2 ICL23 

ai3 —i<223 &33 



(24) 



^(n 5 (5/2)E 

2=1 

+n 5 (3/2)|p? 2 /2 | 2 ) = 



5/2 ,2 



1 



95 -49 cos(2V3Dr) 



J±*(t) = ^n 5 (5/2)|^ 2 | 2 

144 o V7Dt 
■ sin — - — . 



2401 



(29) 



490 L 

- 45 cos(2v / 7^r) - cos(4V7l>t)1 ; 



(30) 
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One can check that the sum of all intensities of Eqs. 
(|29|) , ([30]) equals 1 independently of r according to Ref. 



IV. NUMERICAL CALCULATIONS OF MQ 
NMR DYNAMICS OF A SPIN SYSTEM IN A 
NANOPORE 

The intensities of MQ NMR coherences are oscillating 
functions even at very long times. As an example of such 
behavior one can see the dependencies of intensities Jo, 
J 2l J20, and J22 of MQ coherences on the dimensionless 
time t = Dr which are shown in Fig. [3] for the system 
consisting of N = 201 spins in a nanopore. The mini- 
mal frequency of these oscillations is of the order D and 
the minimal positive eigenvalue of the Hamiltonian Hmq 



obtained from numerical calculations is |A 



(min) I 



V3 



D. 



\3/2) I - 2 

Usually [l| MQ NMR dynamics is closely connected to 
the growth of many-spin clusters which are responsible 
for the emergence of MQ coherences of high orders. In 
the beginning, the clusters consist of spins coupled by 
the strongest DDI which are interactions between nearest 
neighbors. Then, as time goes by, the clusters are grow- 
ing owing to the smaller DDI corresponding to the next 
nearest neighbors. The further growth of the clusters is 
determined by other remote spins. Thus MQ NMR co- 
herences of high orders emerge only at long times. Here 
we have another case. All spins in the nanocavity are 
coupled through the single coupling constant D and the 
cluster consisting of all N spins emerges at the time 
t « 1 (r « 1/D). In fact the process of the growth of 
many-spin clusters is absent in this system. The quasi- 
stationary intensity distribution develops very quickly, as 
it may be seen in Fig. [4j where the distributions taken 
at time moments t = 1,35 almost coincide with each 
other. Significant discrepancies in the beginning of these 
graphs appear due to the oscillating behavior of the inten- 
sities. Nevertheless it takes some time for the emergence 
of the spin clusters which can absorb (emit) n-quanta 
(n = 0, ±2 ± 4, ±6) of the energy of the external rf-fields. 
At small times, an appearance of the corresponding terms 
in the density matrix p can be clarified by the following 
expansion 



p M = e -i{HM Q /D)t Ize i{H MQ /D)t 

I z -itA-^B+ Z -t 3 C + 



(31) 



where 



A = 



(/+) 2 - {I-f 



B = 2/ 2 -I z - 2I z ri+, 



(32) 



(33) 



C = 



-{3(J-) 2 +9(J+) 2 - 
+ 2(/-) 3 /+ 



2/ z (/-) 2 - 14J Z (/+) 2 
-2J-(/ + ) 3 -4J 2 (J-) 2 (34) 
+4I 2 (/+) 2 }/2. 



Below we consider the intensities of MQ NMR coher- 
ences at t > to = 31 when the coherences of all possible 
orders have appeared and one can think that the quasi- 
stationary distribution of the intensities is realized. We 
consider the averaged intensities Jk (-N < k < N) of 
MQ NMR coherences with the averaging over several pe- 
riods T = 2ttD/ |A^ min) i- 



x 3/2 



t +K T 



(35) 



to 



where Kq is some positive integer taken from the require- 
ment that the increase of the averaging interval does not 
change J^. We have found that it is enough to take 
K = 2. 

The distributions of the averaged intensities for TV = 
201, 401 and 601 are shown in FigO One can see that all 
averaged intensities may be separated into two families: 



Ti = {J4k-2,k = 1,2, . . .}, 
T 2 = {J±k,k = 1,2, . . .}, 



(36) 



while Jo may not be reffered to any of them. Each family 
may be approximated by the smooth function so that the 
distribution of the intensities J2n can be represented as 
follows: 

Al(1 + 2ai|n| + 4a 2 n 2 ) exp(-2ai |n|), 
i2n= \ n = ±l,±2,... . (37) 

A 2 exp(-2a 2 |n|), n = ±2, ±4, . . . 

Parameters oli and ai (i = 1,2) have been found for 
N ■ 



201, 401 and 601: 
N = 201 : A x 



(38) 



0.0875, A 2 = 0.0912, 
ai = 0,0838, a 2 = 0.0570, 
ai = -0.0648, a 2 = 0.0059, 
N = 401 : Al = 0.0560, A 2 = 0.0704, 
ai = 0,0577, a 2 = 0.0429, 
ai = -0.0433, a 2 = 0.0030, 
N = 601 : A x = 0.0437, A 2 = 0.0608, 
ax = 0,0465, a 2 = 0.0367, 
ax = -0.0345, a 2 = 0.0020. 

These parameters allow us to fulfill the normalization 
condition 11811 : 



Jo + 2 ^ Jk = 1, 



(39) 



k=l 



where we take the limit N — » 00. 
Eq. (|39|) can be rewritten as 



Jo 



2A 2 



{sinh 2 (2ai) 



e 4« 2 _ 1 sin h 3 (2ai) 
ax sinh(4ai) + 6a 2 + 2a 2 cosh(4ai)} = 1, (40) 




Time,t 



FIG. 3: Oscillating dependences of the intensities of MQ NMR coherences on the dimensionless time t in the system of N = 201 
spins in a nanopore; a) Jo and J2; (b) J20 and Jii- 
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FIG. 4: Evolution of the MQ NMR coherence profile in the 
system of N = 201 spins in a nanopore at t = 0.1, 1, 35. 



where Jo must be found for each particular number of 
spins N. For instance, we have obtained the following 
three values of Jq: 



0.1973, TV = 201 
Jo = { 0.1604, N = 401 . 
0.1414, N = 601 



(41) 



The values of the parameters in Eq. (|37|) have been found 
according to the procedure consisting of the following 
three steps. 

1. The family of the intensities T 2 of MQ NMR 



FIG. 5: The averaged intensities of MQ NMR coherences J n 
over the time interval to < t < to + 2T with to = 31 and 
T — 4tt/V3 in the systems of N = 201, 401 and 601 spins in 
a nanopore. 



coherences is approximated by the function 
A2 exp(— 2^2 M) with two arbitrary parameters A2 
and OLi- 

2. The value of A2 is expressed via all other parame- 
ters using Eq. (|40|) . 

3. The family of the intensities T ± of MQ NMR co- 
herences is approximated by the function A±(l + 
2ai|n| + 4<22n 2 ) exp(— 2ai|n|) with three arbitrary 
parameters ai, a2, and a\. 
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The results of these approximations are represented in 
Fig. [5j We emphasize that the intensities of MQ co- 
herences have exponential dependence on their orders. 
This model does not confirm the Gaussian character of 
the MQ NMR coherence profile [1]. On the contrary our 
analysis yields a similar exponential profile of MQ NMR 
coherences to what was found for MQ NMR coherences 
in solids in Ref. [4|. In addition, the possibilities of the 
suggested method allow us to find also some subtle pe- 
culiarities of the MQ NMR profile which are typical for 
the spin clusters in nanopores. 



CONCLUSIONS 



Presidium of RAS No. 27 "Foundations of fundamental 
investigations of nanotechnologies and nanomaterials" . 
APPENDIX 

In order to prove the equality ([T6]) we use the following 
relationships: 



N/2 

c n+i 

s=o 



N/2 



2 Ar ;£^- 1 = 2JV - 1 -k 



S=l 



N/2 

2 ~n ; 



N/2 

£c*-=2^ 2 -icr, (A.l) 



S=2 



We present the model in which MQ NMR dynamics 
can be analyzed with numerical methods for systems 
consisting of several hundreds of spins. It is essential 
that this model can be experimentally investigated by 
the methods of MQ NMR spectroscopy. Since the av- 
eraged coupling constant D depends on the volume and 
shape of the nanopore and its orientation relatively to 
the external magnetic field [13, EH , one can obtain this 
information from the comparison of the calculated MQ 
NMR spectra with the experimental ones. The developed 
method can be also used in order to investigate MQ NMR 
dynamics in fluctuating nano-containers [19j]. Then the 
information about correlation times of these fluctuations 
can be obtained [19]. 

Usually [1], MQ NMR dynamics is studied on the ba- 
sis of the averaged non-secular two-spin/two-quantum 
Hamiltonian of Eq. [2j At the same time, this approach 
allows us to take into account the corrections due to the 
higher order terms of the averaged Hamiltonian theory 
[201 ], since these terms commute with I 2 as well. MQ 
NMR dynamics can be also investigated for systems with 
dipolar ordered initial state |2l], [22[ by the suggested 
method. 

Our approach demonstrates that, on the one hand, the 
MQ NMR spectroscopy is very useful tool for the investi- 
gations of the nano-size systems and, on the other hand, 
the study of such systems allows us to clarify some subtle 
peculiarities of the MQ NMR dynamics. 



where C# = (£) - SKN _ S)] 



tions yield 



N/2 



N\(2S+1) 



-(^ + 5 + 1)1(^-5)! 



. The simple transforma- 



(A.2) 



N/2 



N- 



N/2 

2 n N/2-S 4 V- q r N/2-S 

N+l jLi D ' °iV+l 



s=o 



s=o 

N/2 



1 r N/2-b 
N+l 2^ °AT+1 • 



One can find directly with Eq.(Al) that 



4 N.W 



__2 N+1 +2C N N I\ (A.3) 



N + l ^ N+1 

s=o 



N/2 



N- 



s=o 



Using Eqs. (|A.1[) . (|A.3|) . (|A.4[) one obtains 



(N-2)-2 N -2C% /2 , 



(A.4) 
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iN 
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The equality of Eq. ([T6]) is proved. 
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